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Abstract 

Even if never elucidated, the question of determinism is a standing question along the 
history of human thinking. A physical system evolves in a deterministic way if its future 
is completely determined once we have fixed some present characteristics of it, i.e., its ini- 
tial conditions. The problem addressed in the present paper is to test determinism in the 
macroscopic domain. By imposing a very plausible "separability" assumption, we prove that 
determinism enters in contradiction with the recorded outcomes of a given electrocardiogram. 
The interest of this result comes from the fact such a basic idea as determinism has never been 
experimentally tested up to now in the macroscopic domain, and as far as we know not even 
in the quantum domain. 
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1 Introduction 

It is well known that realism, the previous existence of some hidden variable values determining 
any given outcome, has been experimentally tested in the quantum domain [T], [2J, [3J. The 
experimental result is that, according to quantum mechanics predictions [3], [S], what is called 
local realism has to be discarded, although some loopholes remain still open [5]. Here, "local" 
means that the corresponding reality -the hypothetical so-called hidden variables- cannot involve 
correlations among its distant enough parts, i.e., among the parts unable to interchange signals 
because of the relativistic limit of interaction propagation, i.e., because of the limit of the speed of 
light. 

Nevertheless, what we will test in the present paper is not realism, but determinism. In order 
to catch the difference between both concepts, consider two consecutive self-outcomes of a given 
isolated system. Imagine that we have some reality behind any outcome (realism) . Let us consider 
the reality which is behind the first outcome. Then, it could happen that the same outcome erases 
this reality such that the reality which is supposed to be behind the second outcome is a different 
one (this is just what happens with the quantum state collapse in the measurement process). In 
such a case, we would have realism, but not necessarily determinism: in order to have determinism 
some initial conditions must exist, the same for all the system evolution. This is what is assumed 
to happen, implicitly at least, in classical physics. 

But, letting aside the interest of determinism as such, what could be the remaining interest 
of testing determinism when realism has already been tested with the resulting rejection of local 
realism? There are two reasons for this interest: in the first place, in the present paper we 
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test macroscopic determinism by dealing with a very ordinary macroscopic system, and such a 
macroscopic test has never been performed as long as we know; in the second place, in the present 
paper framework, the locality condition is a too hard constraint since, even accepting the speed of 
light limit, as it should be, we could have some of the above precluded non local correlations by 
simply allowing the system to evolve in an appropriate deterministic way (a possibility to which 
Bell |7] referred as a "mind boggling . . . conspiratorially entangled" world). Thus, what is really 
at stake behind the question of realism is just determinism and this is why we directly consider it 
here. 

2 The Case of an Electrocardiogram 

More precisely, we have considered the record of the successive times corresponding to the so 
called R self-responses (the typical sharp peaks) of an electrocardiogram. The file, together with 
the code used in the analysis that follows it (see the Appendix) contains 4340 recording times, 
from an anonymous patient's night hour record, whose mean time interval between two successive 
self-responses is tj\/ = 829 milliseconds. Then, we consider the time succession {t n = 3ntM}, 
with n — 0, 1, 2, 3, . . . , 1446. Let us consider too three time intervals r a ,a = 1,2,3, such that 
T i < T 2 < T 3, with t 3 < 3tnj. From the file of the above recording times, we will define a 
dichotomic (i.e., two-valued) function U(t n + r a ) = ±1 in the following way: If r a < tM, U takes 
the value +1 if during the time interval (t n ,t n + r a ) at least a response R appears, and the value 
-1 if it does not appear during this interval. Similarly in the two remaining cases: tu _ T a < 2tM 
with the interval (t n + tM,tn + T a ) 1 and 2£m < t q < Mm with the interval (t n + 2tM, t n + To). 

Then, imagine that for each time t n we choose at random two of the three obtained values 
U(t n + r a ). Let us consider the three corresponding correlation functions 

< U(a)U(b) >= -Tj-Y\U{t n{ab) +T a )U{t n{ab) + Tb ), a<b; a, b =1,2, 3, (1) 

±'ab 

where {t n ( a f,)} is the sub ensemble of the time succession {t n } corresponding to the random selected 
pair (r a , r b ). The summation is over the N ab times i ra ( a f>) and the addition of the three very similar 
numbers N a b will give N . 

In order to test it, let us assume determinism for the appearance of the successive R responses 
and so determinism for the physical system which is behind them. This determinism will be 
assumed, not for the heart producing the electrocardiogram, which is obviously a non isolated 
system, but for some "enlarged system" (the heart plus its actual environment: the entire human 
body, and even the patient's room and beyond) which can be considered isolated o nearly isolated. 
This means, in particular, that some initial conditions, A , referred to some initial time, t , exist, 
such that the above function U can be written more explicitly U(t n r ab -\ + r a ,Ao). The initial 
conditions Ao will collect the values of some variables of the "enlarged system" at this initial 
time. Notice that we also could write U(r a , A n ( ab )) with A„( ;,) denoting the new initial conditions 
corresponding to the times t n ( a b) when these times are taken as new initial times. Obviously, given 
the original initial condition Ao and the new initial time t ra ( a fe), there is a unique initial condition 
A n ( ab ), but this uniqueness is not necessarily valid the other way around: the same value of A 
could correspond to different initial times t n ( ab j with the same Ao. Then, the correlation functions 
of (Q} can be written as 

< U(a)U(b) >= J2Pab( A n(ab))U(T a ,A n{ab) )U(T b ,A n{ab) ) (2) 

where the summation is now over all different A„( a f,) and we have introduced the normalized 
probabilities Pab{A n [ ab )) since the different values of A„( a b) will not appear, in principle, with the 
same frequency. 
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3 Proving some Time Bell-like Inequalities 



Then, we make the following "statistical separability assumption" : the history of "enlarged system" 
(for example, the times in which the R responses appear), is statistically independent of the times 
t n (ab) selected by the authors. This means that in © we can drop the index ab affecting the 
probabilities p a b, and the index n(ab) affecting the A values {A„( ab )}. Thus, we can write as 

< U(a)U(b) >=J2 P(A)(7(r Q , A)U(r b , A) (3) 

Dropping the last index, n(ab), needs some further comment. Notice first that we could write 
(j3|) if the three sets {A„( ab )}, irrespective of the order its elements, were the same set, let us say 
{A}. Trying to ensure this, we could assume that the history of our "enlarged system" is completely 
independent of the ulterior choose of the times £„( a (>) (of course, though stronger, this plausible 
assumption does not enter in contradiction with the above statistical independence) . The problems 
is that, even with this full independence, it could happen that the three sets {A n (ab)} were not the 
same set: i.e., given an element of one of these three sets, it could happen that this element did 
not belong to some one of the other two sets. This could happen if the elements of {A n( - ab )} did 
never repeat themselves. Nevertheless, the function U(r a , A- n (ab)) is a two valued function. This 
means that its values repeat largely themselves. Thus, from the full independence we can expect 
that given any function value, U(r a , A n ( &)), there always will exist a value {A„( ac )} such that 
U(r a , A„( oc )) = J/(to, A n ( a i,)), that is, such that we can write both as t/(r , A) with A belonging to 
the above common set {A}, all which leads to ([3]). However, in order to write © we do not need 
to make such a natural assumption: it is enough to merely assume directly, as we have done, the 
less restrictive condition ([3]) itself because of the assumed statistical independence. 

From the two assumptions, determinism and statistical separability, we next prove the following 
"time Bell-like inequalities" : 

D = | < U(1)U(2) > - < 17(1)17(3) > |+ < U(2)U(3) >< 1, (4) 

where the vertical bars mean taking absolute values. In order to prove inequalities ((4]) we will 
mimic here the proof of the original Bell's inequalities [4j. Having in mind ([3]), i.e., having in mind 
the assumed statistical independence, let us consider the difference 

< 17(1)17(2) > - < 17(1)17(3) >= 5>(A)*7(Ti,A)[/(T 2 ,A) - ^p(A)[/(n, A)U(r 3 , A). (5) 

Notice that a similar statistical independence is implicitly assumed when proving the original 
Bell's inequalities [4] and further similar inequalities: there, the hidden variables, A, are assumed to 
have a random distribution independent of the successive chose of the different pairs of measurement 
directions. A similar statistical independence is also assumed in 8], where it is called the "induction 
postulate" . 

Having in mind that U(r a , A) 2 — +1, inequality ([5]) becomes 

< 17(1)27(2) > - < U(1)U(3) >=Y,P(A-)U(ti,A)U(t 2 ,A)[1-U(t 2 ,A)U(t 3 ,A)}, (6) 

Then, taking absolute values, 

< 17(1)17(2) > - < U(1)U(3) > | < J>(A)[1 - U(t 2 ,A)U(t 3 ,A)], (7) 

that, according to J^p(A) = 1, becomes inequality (|4]), as we wanted to prove. 

Similar inequalities have been proved in the literature on the subject under several assumptions 
different from determinism. One of these assumptions (the non- invasive measurability) [8| [9] being 
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very dubious and the other ones (joint realism |10j and counterfactual macroscopic definiteness [8]) 
unnecessarily strong. In it is claimed that inequalities (|4|) have been proved, in the microscopic 
domain, from the assumption of determinism plus non contextuality. Something similar is claimed 
in |12) . without explicitly requiring the separability assumption. But, as we have just shown, the 
right proof of inequality Q needs this last assumption, besides determinism. 

4 Testing the Violation of the Time Bell-like Inequalities. 
Conclusions 

The main result of the present paper is the violation of inequality (j4]) by the recorded data of 
a given electrocardiogram. This entails the failure of determinism as far as the evolution of the 
above "enlarged system" , including the corresponding heart, is concerned, and so the failure of 
determinism in the natural world. Certainly we could elude this conclusion, but only at the high 
price of allowing the existence of some convenient correlations between the recorded data and our 
arbitrarily differed choice of the times t n and r a . 

The correlation functions < U(a)U(b) > have been estimated using ([T]) for different combi- 
nations of r a ,Tb, a < b, (a,b) — 1,2,3. These are the 6 3 = 216 combinations resulting when 
ti is chosen among the six times (300, 400, . . . , 800), and similarly for T2 among (tM + 300, tm + 
400, ...t M + 800) and t 3 among (2t M + 300, 2t M + 400, ...,2t M + 800). For 19 of these 216 combi- 
nations inequality (0| is violated. Its greatest violation corresponds to the combination t\ — 800, 
7"2 = tn + 600 and T3 = 2tu + 300 for which the value of D in inequality ((4]) becomes D = 1, 12. 

It is necessary to remember that this violation value is produced by a given random choice 
of the sequence of pair times (r a ,rfc), a < b, (a,b) — 1,2,3, for each element of {t n = 3ntM}, 
with n = 0, 1, 2, 3, ... , 1446. There are an enormous number, 3 1446 , of such possible choices. This 
means that the probability of getting the above result, jointly with the other 18 violations, by mere 
chance is very small, as we discuss below. To this purpose, we have generated 10 4 different random 
sequences for the same combination of r's: t\ = 800, = tM + 600 and T3 = 2tu + 300. A total 
of 2251 violations, ranging from 1,01 to 1,24, were encountered confirming that the original result 
is not due to chance in the particular selection of one of the above 3 1446 different sequences. A 
Kolmogorov-Smirnov test has been performed for testing the normality of these 10 4 observations. 
The statistic of the test is 0,0042 with a p-value of 0,4995. So the normality of the data can be 
accepted (Fig. [T]). 

We now confirm, with a deeper analysis, that the above violation is not due to a particular choice 
of the above combination of r a , by significantly shifting this combination t\ = 800, T2 = tM + 600 
and T3 = 2tM + 300, in its neighbourhood. Each time, we have replaced some of the r a by 
r a ± 20. The resulting new combinations have been explored as described above, i.e., by randomly 
generating each time 10 random sequences. This complementary study confirms the original 
violation of D < 1 . In particular, for the combination t\ — 780, T2 = tM + 620 and T3 = 2±m + 320, 
we have even found a greater violation than before, corresponding to D = 1, 28. 

5 Final Considerations 

As we have just explained, this violation entails that the R responses of our electrocardiogram 
have not been generated in a deterministic way, unless we are disposed to leave an assumption 
as unavoidable as the statistical separability. Notice that this assumption does not invalidate 
determinism: our choice of the times t n and r a could be perfectly deterministic, if we want. All 
is needed here is that this arbitrarily delayed choice be independent, or which is less restrictive, 
statistically independent, of the previously recorded R responses. 
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Figure 1: Histogram and adjusted Normal density for the 10 4 repetitions of t's combination giving 
maximum D. 



Notice that we have not needed to assume locality: our assumed initial conditions are not 
required to be local. Assuming locality would have been no consistent in the present case, since 
as we have pointed above one could have non local realism by allowing the system to evolve in a 
deterministic way, and it is just determinism what we have tentatively assumed here in order to 
test it. At the same time, the inefficient detection loophole [3] is not present here, since all the 
self-responses, R, of our electrocardiogram have been considered in our analysis. 

Let us make another remark: the expression for < U(a)U(b) > given in says that we 
deal with a common probability space, i.e., we use the same probability function, p(A), and we 
can put A instead of {^ n (ab)}: all along the history of the system. This common space exists 
because of the particular statistical independence entailed by the separability assumption, which 
has allowed us to drop the index ab and n(ab) from the probabilities p a b and {A„( {,)}, respectively, 
in ([2]). However, could not it be the case that successive new initial conditions were incorporated 
to the actual initial conditions of the R responses along the history of the heart and its assumed 
environment (the "enlarged system"), these new initial conditions not being before present because 
of the relativistic limit of propagation of interactions? If this were the case, it would not be evident 
that such a common probability space exists. But, in the present case, there is not such a problem 
since, because of determinism, there is no exterior effects on this enlarged system, containing the 
heart and its self-responses (a system enlarged as to cover the entire universe if needed), neither 
there is really anything new happening in this system since everything what happens there comes 
from the initial conditions. As a result we must have a common probability space for the A values 
in U(r a ,A) and p(A), such that we be able to write the correlation functions < U(a)U(b) > in 
the way we have done in (|3]). The question of the possible lack of a common probability space 
has also been raised [T3] in the similar scenario of the seminal Bell's inequalities [1], where, as we 
have remarked above, assuming the existence of this space is actually an implicit, even if plausible, 
hypothesis that we do not need to make here because of the determinism assumption. 

Finally, we do not have at present any theory leading to the failure of determinism, out of 



5 



quantum mechanics. Thus, we could speculate that the reported violation of macroscopic deter- 
minism could be the macroscopic magnification of the quantum indeterminism of the microscopic 
constituents of the heart and its environment. Furthermore, we could ask ourselves if this hypo- 
thetical magnification comes from the fact that we were dealing with living stuff, or more generally 
it would be expected to happen every time that we have some sort of dichotomist self-response in 
an, living or not, isolated system. All these questions would deserve some future work. 
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Appendix 

Readers can recover the data file and can reproduce the analysis results, described in the article, 
downloading the files ecgn.txt and file analysisR.txt from http: / /www. uv.es/montes/electrocardiogram 
. The first file contains the data corresponding to the electrocardiogram with individual and ac- 
cumulative response times, the second one is an R [14j code allowing to reproduce the analysis 
performed. This code provides three output files: 

outputl.tx: the results generated by the 216 combinations of t±, t 2 and T3. The columns are, in 
this order, n, r 2 - t M , t 3 - 2t M , D, ni 2 , n\z and n 23 

output2.txt: the results generated by 10 4 random repetitions obtained for the combination of t's 
providing the maximum D among the 216 of the above file. The file has the same structure 
as outputl.txt and, obviously, n, t 2 and T3 have the same value along the 10 4 repetitions. 

summary.txt: a self explicative file summarizing the analysis, including also the result of the 
Kolmogorov-Smirnov test for testing the normality of the 10 4 D values in the above file. 

It is possible to obtain an exact copy of our results if the same random seed is used at the 
beginning of the R code. In our case this value was 12345. 
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